Controllable plasma energy bands in a ID crystal of fractional Josephson vortices 
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We consider a ID chain of fractional vortices in a long Josephson junction with alternating i^c 
phase discontinuities. Since each vortex has its own eigenfrequency, the inter- vortex coupling results 
in eigenmode splitting and in the formation of an oscillatory energy band for plasma waves. The 
band structure can be controlled at the design time by choosing the distance between vortices or 
during experiment by varying the topological charge of vortices or the bias current. Thus one can 
construct an artificial vortex crystal with controllable energy bands for plasmons. 
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The study of crystals is a cornerstone of solid state 
physics 0,13,0. The electronic structure of crystals, such 
as energy bands and the dispersion relation for electrons 
and phonons, has a great impact on all physical proper- 
ties of solids. A crystal is a periodic arrangement of a 
group of atoms translated in space by direct lattice vec- 
tors (usually edge effects are negligible). This periodicity 
in space results in the formation of energy bands and dif- 
ferent propagation modes for electrons and phonons. If 
one could control the periodicity, e.g. by varying the mu- 
tual position and the kinds of atoms in a crystal, during 
experiment, one could turn a metal into an insulator and 
basically create a material with any desired properties. 
Unfortunately, the crystal structure (the type of ions, lat- 
tice and spacing) is fixed by nature and can be changed 
only a little, e.g. by applying strong pressure, electric or 
magnetic fields. 

Therefore, it is interesting to construct artificial peri- 
odic structures (crystals), whose properties can be varied 
over a wide range during experiment^ allowing a great 
degree of control over their resulting electronic prop- 
erties. Such structures must not necessarily consist of 
atoms, but could consist of other more macroscopic ob- 
jects which can be fabricated and arranged by means of 
modern technology, e.g.^ using lithography. 

In this letter we consider an artificial ID crystal made 
of fractional Josephson vortices (fractional magnetic flux 
quanta), and study its plasma energy bands, ie., energy 
bands of plasmons — small oscillations of the Josephson 
phase. 

Usually, vortices in Josephson junctions and supercon- 
ductors (S) carry one quantum of magnetic flux <I>o ~ 
2.07 X 10~^^Wb. Fractional Josephson vortices spon- 
taneously appear in long Josephson junctions (LJJs) 
with an unusual current-phase relation 0, I^, with spa- 
tially alternating sign of the critical current {e.g., 
SFS LJJs0, 1^ with a step- wise variation of the fer- 
romagnetic (F) layer thickness) or with tt discontinu- 



ities of the Josephson phase i^| (e-q., dx2_y2-wdive based 
LJJs[i3, m El El El El [S Ell). More general, a 
/^^-discontinuity of the phase at x = (x is the coordi- 
nate along the LJJ) means that the relation between the 
supercurrent Ig and the phase ja (the first Josephson re- 
lation) in the region x < is Ig = /csin(/i), while in 
the region x > it is = /c sin(/i + k). Fractional vor- 
tices spontaneously appear to compensate the ^-phase 
jump and are pinned at the discontinuity point sjl^- For 
a TT discontinuity the vortex carries the fractional flux 
^0/2 and the junction can be described by the same 
model as the junction with alternating sign of critical 
current HH, El El m El mentioned above. 

Recently, a LJJ with artiflcial phase discontinuities was 
proposed, implemented and successfully tested[2l. In 
this junction, made using standard Nb-AlO^^-Nb tech- 
nology, the discontinuity of the phase is created using a 
pair of tiny current injectors attached to the same elec- 
trode of the junction as close as possible to each other. 
By passing the current /inj from one injector to the other, 
one can create an arbitrary /^-discontinuity of the Joseph- 
son phase with hi oc /inj. Since the Josephson phase is 
27r-periodic, without loss of generality we consider only 
< < 27r. One can also fabricate such a junction 
with as many injector pairs as required and tune the 
strength of all discontinuities by a single control cur- 
rent. Nb-A10a;-Nb junctions also have very low damping 
(McCumber- Stewart parameter /3c ^ 1), which makes 
them a perfect candidate to study propagation of plasma 
waves and formation of energy bands. In comparison, 
SFS and d-wave based LJJs typically have /3c ^ 9 (see, 
e.g., Fig. 2 of Ref.lTil). 

A fractional vortex has an eigenfrequency c<;o(/^) which 
corresponds to oscillations of the magnetic flux around 
the discontinuity point ^2^. This eigenfrequency lays 
within the plasma gap and depends on the flux carried by 
the vortex [2^. If one considers two vortices at some dis- 
tance from each other, so that their magnetic flelds and 
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super currents overlap, then the vortices behave as two 
coupled oscillators: their eigenfrequency splits into two 
frequencies, corresponding to two different modes: in- 
phase and out-of-phase. This splitting of eigenfrequen- 
cies for different two-vortex-configurations was already 
investigated numerically [2^. In an infinite ID array of 
fractional vortices situated not very far from each other, 
due to coupling, the eigenfrequencies split and form an 
energy band, very similar to the formation of bands in 
a crystal starting from the discrete energy levels of elec- 
trons in a single atom. In this paper we study the energy 
bands corresponding to oscillations of a chain of frac- 
tional vortices and we show that the band structure can 
be easily controlled during experiment. 

Consider an infinite array of alternating disconti- 
nuities at a distance a from each other and ^k, vortices 
pinned at them. Such an antiferromagnetically (AFM) 
ordered state represents the most natural configuration 
or, in other words, the ground state of the system. 

To model an infinite chain of vortices we consider an 
annular LJJ with two phase discontinuities placed at the 
same distance a from each other, so that 2a = L, where L 
is the length (circumference) of the LJJ. It is also possible 
to realize such a geometry experimentally using an annu- 
lar LJJ with injectorsjl^, I23, E^- The dynamics of the 
Josephson phase is described by the following perturbed 
sine- Gordon equation 0, 0, 

- Mtt = sin [/i + e{x)] - 7, (1) 

with 

eix) = l 0' 0<x<a, 

^ ^ I -K, a < X <2a = L. ^ ^ 

The Josephson phase ii{x^t) is a continuous 2 7J func- 
tion of the spatial variable x and time t, normalized to 
the Josephson penetration depth Aj and to the inverse 
plasma frequency oo~^^ respectively. The subscripts x 
and t denote derivatives with respect to space and time. 
We normalize the bias current density j to the critical 
current density jc, i.e., 7 = j/jc The function 9{x) 
defines the positions and the strength of the discontinu- 
ities. Since we consider an annular geometry we sup- 
plement Eq. (QJ with the periodic boundary conditions: 
/i(0,t) = and /^x{0,t) = iix{L,t). 

To calculate numerically the oscillatory energy bands 
of a chain of fractional vortices, we, first of all, cal- 
culate the steady solution corresponding to an AFM 
vortex chain for given values of a, and 7 by solv- 
ing Eq. Q numerically ^2^. We use several spatial dis- 
cretization steps to compare the obtained results, i.e.. 
Ax = 0.1, 0.05, 0.02. 

If /io(^) is the obtained static solution of Eq. (QJ, then 
the small oscillations around this solution are given by 

IJL{x,t) = /j.o{x) ^ e{x)e''^K (3) 



Substituting expression ^ in Eq. we get 

+ Lj'^e = cos[/io(^) + 0{x)]e. (4) 

Written as a system of first order equations, the eigen- 
value problem takes the form 

Cx = Ac, where ^= ^ ^ ^ (5) 
with the matrix 

^^^^ " (cos [/io(x) + 0{x)] -00^ 0) • 

In our case, since A{x) 7^ const^ we cannot implicitly 
integrate Eq. (0). Instead, we approximate A{x) by a 
constant matrix An given by Eq. Q with x = Xn = 
nAx on each interval Ax from Xn to Xn-\-i- Then we 
solve Eq. © on each small interval Ax to find e{xn-\-i) = 
e{xn)exp{AxAn). 

As a consequence of the Floquet (Bloch) theorem the 
solution of Eq. Q should satisfy the following equation 

e{x + 2a) = Ce(x), 

with constant matrix C, where 2a is the periodicity of the 
lattice. In our case, the so-called principal matrix C is 
approximated by a simple transfer matrix discretization 
described above, i.e., 

N 

C= nexp(AxAn), (7) 

n=l 

where N = 2a/ Ax. More delicate schemes to calculate 
a band-gap can be found in, e.^., Refs. jHH^. The 
transfer matrix C has two eigenvalues. The product of 
these eigenvalues 

N N 

A1A2 = detC = Y[ det e^^^- = Y[ e^'^^^^-^ = 1, 

n=l n=l 

where we used an identity det[exp(M)] = exp[Tr(M)] 
and the fact that Tt{A) = 0, see Eq. The periodic 
solution exists, i.e., uj is within the energy band, only if 
I All = IA2I = 1, i.e., if the determinant of the character- 
istic equation is negative and Ai^2 is a pair of complex- 
conjugate roots laying at the unit circle on the complex 
plane. 

First, we consider the case when the bias current is 
absent, i.e., 7 = 0. Band structures as a function of 
K,^ numerically calculated for different a (given in units 
normalized to Aj) from a = 5 (weak coupling) to a = 1 
(strong coupling) are shown in Fig. ^ In all plots one 
can see that in the absence of discontinuities {hz = 0) 
the junction has a plasma gap for < < 1 (cj is nor- 
malized to the plasma frequency cOp) and a single infinite 
plasma band for a; > 1. As increases, fractional vortices 
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FIG. 1: Numerically calculated band structure as a function of n for a = 5, 4, 3, 2.5, 2, 1 shown in (a)-(f), respectively. At 
n = Kc(a) a crystal becomes unstable and structural phase transition takes place. 



appear. Each vortex, if isolated, has an eigenfrequency 
^o(^) < 1^3- our case vortices are coupled and the 
eigenfrequency splits into a band, which is the lowest 
energy band, see Fig. ^a). Small gaps also appear in 
the former continuous plasma band a; > 1. As the cou- 
pling increases (distance a decreases), the bands broaden, 
while the gaps shrink and shift to higher frequencies, as 
can be seen in consecutive Figs, n^a)-(f). 

In Fig. n the bands are traced from = to = 
Kc{a). At = nc{ci) the AFM state becomes unstable 
and turns itself into a complimentary state, i.e., each of 
the vortices becomes a — 27rVvortex|l8j. As a 
sign of this instability the lower band touches zero, ie., 
uo ^ {) K ^ Kc- For the complimentary state, the 
band structure is the mirror reflection of the one shown 
in Fig. Q with respect to the line = ir. Note, that in the 
interval of hi from 27t — hcda) to K^do^) there are two stable 
solutions: an AFM chain of direct vortices or an AFM 
chain of complementary vortices. Similar behavior was 
reported for a system of only two vortices 0. The value 
of K,c{o^) decreases as the coupling increases reaching the 
value ^ TT at a ^ 0, in agreement with the known 



result that the infinite AFM ordered semifluxon 
chain is stable for any a ^ 0. 

Contrary to a crystal of integer vortices j^E 1^ 135^ . 
our ID crystal has no acoustic branch (except the case 
n = 27rn) in the dispersion relation because a^- vortices are 
pinned and cannot move all together, i.e., with 00 = and 
k = 0. 

Next, we consider the control of the band structure by 
the bias current at various fixed values of a and = ir. 
This case is interesting for O-tt-L J Js where the discontinu- 
ities cannot be controlled during experiment. In Fig. Owe 
present the band structure for the case of weak (a = 5), 
moderate (a = 2) and strong (a = 1) coupling. As one 
can see, by applying a bias current, an additional gap 
opens within each band. For the case of moderate and 
strong coupling these gaps can become quite large near 
the critical value of 7 = 7c (<^) at which the static solu- 
tion becomes unstable. Again, as a sign of instability the 
lowest band has cc; ^ at 7 ^ 7c(a). For 7 > 7c(«) the 
system switches to the finite voltage state. Note, that 
7c (a) ^ at a ^ in qualitative agreement with the 
results for finite semifluxon chains. g31i] 
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FIG. 2: Numerically calculated band structure as a function of the applied bias current 7 for At = tt and (a) weak coupling 
a = 5, (b) moderate coupling a = 2 and (c) strong coupling a = 1. At 7 = 7c (a, At) a crystal becomes unstable and structural 
phase transition takes place. 



To conclude, we have calculated the energy bands cor- 
responding to small oscillations (plasma waves) of the ID 
AFM ordered fractional vortex crystal as a function of the 
discontinuity hz. Such a ID vortex crystal has no acous- 
tic, but only an optical branch in the dispersion relation, 
which is a direct consequence of vortex pinning. The 
band structure can be changed by changing k,. In case 
of discontinuities created artificially by injectors, one can 
make a wiring such that a single control current changes 
the value for all discontinuities in the same time, thus 
providing the possibility to chan ge t he band structure 
"on-the-fly" . For natural O-tt-L J Js [ll E [T^ , where the 
discontinuity = ir is fixed, the band structure can be 
smoothly controlled during experiment by the bias cur- 
rent. 

The knowledge of eigenmodes and the band structure 
is a key element in designing classical or quantum devices 
based on fractional vortices. In the classical domain this 
may help to avoid resonance phenomena or to exploit 
them {e.g. in filters and detectors). In the quantum 
domain, thanks to the absence of acoustic branch, the 
lowest energy gap can be a crucial obstacle for thermal 
excitation of plasmons. 

Here, we have calculated the plasmon spectrum for a 
mirror symmetric crystal. However, the most unusual 
properties can be expected from systems with broken re- 
flection symmetry (ratchets) [sE 5 such as crystals of 
ferroelectrics or of some superconductors. From this 
point of view, the transport in such systems is not well 
studied. Using an array of discontinuities of different 
strengths and distances, one might realize controllable 
fractional vortex crystals without reflection symmetry 
and study the nonequilibrium transport. In this case 
the eigenvalue problem corresponds to the motion of 
a plasmon in a ratchet potential. 
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